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Abstract. We study first-passage time problems for a diffusive particle with 
stochastic resetting with a finite rate r. The optimal search time is compared 
quantitatively with that of an effective equilibrium Langevin process with the same 
stationary distribution. It is shown that the intermittent, nonequilibrium strategy 
with non-vanishing resetting rate is more efficient than the equilibrium dynamics. Our 
results are extended to multiparticle systems where a team of independent searchers, 
initially uniformly distributed with a given density, looks for a single immobile target. 
Both the average and the typical survival probability of the target are smaller in the 
case of nonequilibrium dynamics. 
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1. Introduction 

Stochastic search problems occur in many fields of science as well as in daily life. The 
quest for an optimal strategy for locating a target whether inanimate (such as a binding 
site for a protein at the molecular level pQ or an element in a list) or living (such as a 
prey for a predator [2]) has been the source of a large number of different algorithms 
that combine observation, physical mechanisms and computation. Depending on the 
context, search strategies can be very different, leading to a variety of interesting 
models (see the special issue [3] devoted to this field of research). One robust class of 
models called intermittent target search strategies combine phases of slow motion, that 
allow target detection, and phases of fast motion, during which the searcher relocates 
but is not reactive (see [I] for a recent review). Such strategies have been observed 
at different scales: foraging animals, such as humming birds or bumblebees, display 
intermittent search patterns [HE]. The E. Coli bacteria alternates ballistic moves (or 
'runs') with random changes of direction ('tumbles') in order to reach regions with 
high concentration of a chemo-attractant (chemotactic search) [H EJ |9]. A protein 
efficiently localizes a specific DNA sequence by alternating Id sliding phases with free 3d 
diffusion ('relocation phases'): this mechanism of 'facilitated diffusion', first proposed 
by Adam and Delbriick [1] , enhances the association rate by two orders of magnitude as 
compared to the diffusion limit and leads to numerical values that agree quantitatively 
with experimental results (TQl [TTJ (see [12] for a review). 

A simple model of diffusion with stochastic resetting, in which a Brownian particle 
is stochastically reset to its initial position with a constant rate r was defined and 
studied in [13]. The stationary state of this process is a non-Gaussian distribution and 
violates detailed-balance: a non-vanishing steady-state current is directed towards the 
resetting position. This process can be viewed as an elementary model of an intermittent 
strategy in which the searcher, having explored its environment unsuccessfully for a 
while, returns to its initial position and begins a new search. It was also shown in [13] 
that there exists an optimal resetting rate r* that minimizes the average hitting time to 
the target. Extensions to space depending rate, resetting to an random position with a 
given distribution and to a spatial distribution of the target were considered in [14J. 

The effect of resetting was previously studied in a stochastic multiplicative model 
of population growth where stochastic resetting events of the population size was shown 
to lead to a stationary power-law distributed population size distribution [T5]. A 
continuous-time random walk model in the presence of a drift and resetting has also 
been studied recently [TB]. Finally, in the context of search process, a related model 
has been studied by Gelenbe [T7] where searchers are introduced stochastically into 
the system: there is a single searcher present at a given time with a random lifetime 
and when the searcher dies, a new searcher is introduced into the system at the initial 
starting point. 

In the mathematics literature, mean first-passage time for a class of random walks 
with stochastic restarting events has been studied recently from an algorithmic point of 
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view [T8] . 

In the present work, we focus on the simple continuous-time Brownian diffusion 
model with a nonzero resetting rate introduced in [13] and compare quantitatively the 
reset dynamics with the equilibrium Langevin process that leads to the same stationary 
distribution. We prove analytically that the non-equilibrium process is more efficient 
than the equilibrium one by showing that the optimal mean-first passage time is smaller 
in the former case. We extend our study to the multiparticle problem by considering a 
team of searchers uniformly distributed on the line and looking for an immobile target. 
We show that both the average survival probability and the typical survival probability 
of the non-equilibrium process are smaller than the corresponding quantities for the 
equilibrium dynamics. This shows that, at least for this model for which exact analytical 
calculations can be performed, the non-equilibrium strategy systematically defeats the 
equilibrium behaviour. 

The outline of this work is as follows. In section |2j we review the diffusion model 
with resetting and some of its basic properties, including survival probability and first- 
passage time. An effective equilibrium dynamics with the same stationary distribution is 
defined in section [3] and it is shown that the optimal mean first-passage time associated 
with that dynamics is always larger than the optimal time if resetting is allowed. 
In section [4], we investigate the survival probability of a single target in presence of 
a team of independent searchers. Average and typical survival probabilities display 
drastically different behaviours (algebraic vs exponential decay with time) both for the 
nonequilibrium and equilibrium processes. Nevertheless, it is shown that nonequilibrium 
dynamics is systematically more efficient. The last section is devoted to concluding 
remarks. 

2. Diffusion with stochastic resetting: nonequilibrium dynamics 

In this section, we recall the definition of the model of diffusion with stochastic resetting 
and briefly review some basic results, derived in [13J. 

We consider a single particle on an infinite line starting at the initial position xq at 
t = 0. The position of the particle at time t is updated in a small time interval dt by 
the following stochastic rule [T5] : 

x(t + dt) = xq with prob. r dt 

= x(t) + n(t)dt with prob. (1 - r dt) (1) 

where rj(t) is a Gaussian white noise with mean (i](t)) = and the two point correlator 
(r](t)r](t')) = 2 D 5(t—t'). The dynamics thus consists of a stochastic mixture of resetting 
to the initial position with rate r (long range move) and ordinary diffusion (short range 
move) with diffusion constant D (see Fig. (jl|). Resetting introduces a new length scale 
«q 1 = y/D/r in the ordinary diffusion problem. 

As shown in [13], the probability density p(x,t) of the particle evolves via the 
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resetting rate = r 
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Figure 1. Schematic space-time trajectory of a one dimensional Brownian motion 
that starts at xq and resets stochastically to its initial position xq at rate r. 



following Fokker-Planck equation 

d t p(x, t) = D dlp(x, t) — r p(x, t) + rS(x — xq) 



(2) 



starting from the initial condition p(x, 0) = 5(x — xq). We emphasize that the dynamics 
violates detailed balance manifestly: the current from a position x to Xq via the resetting 
move is not compensated by a current from xo to x. 

The Fokker-Planck equation ^ admits a stationary solution in the t — > oo limit, 
given by 

(3) 



Pst{X) 



Q'o 



exp [— «o | a; — xq\] where ao = yr/D . 



Even though the stationary solution can be expressed as an effective Boltzmann weight: 
Pst(x) oc exp [— Kfrt^)] with the effective potential V e $(x) = a \x — x \, one should keep 
in mind that this is actually a nonequilibrium stationary state, and not an equilibrium 
stationary state. Indeed, in this stationary state detailed balance is violated by a nonzero 
current in the configuration space. Had the stationary state been an equilibrium one, 
the current would have exactly vanished. 

As shown in [131 D3]> t ne additional resetting parameter r allows us to tune the 
target search process to make it more efficient. Consider an immobile target at the 
origin x = and let the searcher undergo diffusion with resetting dynamics specified in 
Eq. ([T]) starting from the initial position xq > 0. What is the mean first-passage time 
T reset (xo, r) to the origin knowing that the searcher starts at x$l This average time, 
T rese t(xo, r), that the searcher takes to find the target will be taken as a measure of the 
efficiency of the search process: the smaller the value of T rese t(^o; r ) (f° r a fixed x ), the 
better is the search. 
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In order to calculate T reset (xo, r), it is useful to consider the survival probability 
q rese t{xo,t) of the target given that the searcher starts at the initial position xq and 
resets to Xq. The mean-hitting time is then given quite generally by [20], ETJ [22] . 
T reset (x ,r) = / °° t (-d t q Teset (x , t)) dt = f °° q reset (xo,t) dt. 

More generally we define Q(x,x ,t), the survival probability of the target given 
that the starting position of the searcher is x and its resetting position is xq. The initial 
position x of the searcher can be treated as a variable leading to a backward Fokker- 
Planck equation for Q(x, xq, t) [131 EH]- Eventually, at the end of the calculation, one sets 
X = Xq and obtains q(xQ,f) = Q{x = XQ,Xo,t). The backward Fokker-Planck equation 
reads for x > 

B t Q(x, Xq, t) = D d%Q(x, Xq, t) - T Q[x, Xq, t) + T Q(x , Xq, t) . (4) 

We require the initial condition Q(x,XQ,t = 0) = 1 for all x > and impose the following 
boundary conditions: Q(x = 0, xq, t) — and Q(x — > oo, xq, t) is finite. Equation |4]) can 
be solved explicitly in terms of the Laplace transform Q(x, x , s) = J °° Q(x, x , t) e~ st dt. 
In particular, setting x = Xq, we obtain 

1 - exp (-J^xoj 

g re set(xo, S) = Q(x = Xq, Xq, s) = ^ \ , (5) 

s + r exp [ — J^jfxo 



thus leading to the mean first-passage time 

f°° 1 

^resetOo, = / "?( X 0, *) dt = q(x , S = 0) = - [exp(« Xq) ~ 1] , (6) 
JO r 



where we recall that ctQ = y/r/D. In terms of the dimensionless parameter 7 = ao^o 
yrjD xq, we obtain 

e 7 - 1 

Preset (^0,7) = 



o/2 



•' 



< 7 ' 



We observe that T rese t(^o, 7) diverges in both the limits 7 — > (pure diffusion) and 
7 — > 00 (for an infinite resetting rate the particle is localized at Xq). The mean first- 
passage time has a unique minimum (see Fig. (J2])) at 7 = 71 where 71 is the solution 
of <9 7 T = 0, i.e., 71 = 2(1 — e~ 71 ), giving 71 = 1.59362 .... Hence, for fixed values of 
D and Xq, there is an optimal resetting rate r = r% = , yfD/x l that makes the search 
time minimum and the search process most efficient. The corresponding optimal mean 
first-passage time is given by 



rpOpt 



reset 



e 71 - 1 



X _ 1 CAA 1 A X 



D 1.54414...-^. (8) 



3. An Effective Equilibrium Dynamics 

We have seen that the resetting with diffusion leads to a current- carrying nonequilibrium 
stationary state given in Eq. ([3]). This stationary distribution can be expressed as a 
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Boltzmann weight p s t{x) = («o/2) exp [— V^x)] with an effective potential 

V eS (x) = a \x - x \ , (9) 

where a was defined in Eq. 

Let us now consider the following Langevin evolution of the particle position with 
time t 

dx 

— = -Bd x V cS (x)+r ] (t) (10) 

where B is the amplitude of the external force and r](t) is the same Gaussian white 
noise as in the previous case, i.e., with {rj(t)) = and (r)(t)r](t')) = 2D5{t — t'). Note 
that we have chosen the same diffusion constant D as in the resetting case to reflect the 
fact that without the reset (in the former case) or without the external potential (in 
the Langevin case), this equation describes ordinary diffusion with the same diffusion 
constant D in both cases. The corresponding Fokker-Planck equation for the probability 
density P(x,t) of the particle is 

d t P(x,t) = Dd 2 x P(x,t) + d x [B (d x V cS (x)) P(x,t)\ . (11) 



Eq. (11) can be expressed as a continuity equation, d t P = —d x J where the current 
density J(x, t) = —Dd x P—B (d x V e g (x)) P. The system will then reach a stationary state 
at long times and if we set the current in the stationary state to be 0, we arrive at the 
equilibrium Gibbs-Boltzmann solution which reads, P eq (x) = Af exp [— (B/D) V e g(x)}, 
where M is the normalization constant. By choosing B = D and V e s(x) = ao \x — xq\, we 
can engineer a zero current equilibrium state which has the same weight as the current 
carrying nonequilibrium stationary state in the resetting case, i.e., 

P eq (x) = p st (x) = y exp [-a \x - x Q \] where a = \JrjD . (12) 

The following natural question then arises. Consider the target search problem, 
where we have an immobile target at the origin. In the previous section, the searcher 
was performing normal diffusion with stochastic resetting to its initial position xq- Now, 



suppose that the searcher undergoes instead the Langevin dynamics as in Eq. (10) 
with the choice B = D and V e g(x) = ao\x — Xq\ that guarantees that both dynamics 
lead to the same steady state. One is then tempted to compare the efficiency of the 
nonequilibrium search process with diffusion and reset to the Langevin search process 



where the searcher's position evolves via Eq. (10). Which process is the more efficient? 



To address this question, we need to compute the mean first-passage time 



Tia, age (xo, r) to the origin of the Langevin process in Eq. (10). With the choice B = D 



and V e s(x) = «o \x — xq\ where «o = \frJP ) 1 the Langevin equation reads 

dx 

— = -a D sgn(x - x ) + rj(t) (13) 
where sgn(^) denotes the sign of z. We define Q\a, n ge(x,x ,t), the probability that the 



searcher, starting at the initial position x > and evolving via Eq. (13), does not 
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reach the origin (the target) up to time t. Equivalently, Q\g, n ge(x, xq, t) is the survival 
probability of the target up to time t under the Langevin dynamics of the searcher. 
Treating the initial position x as a variable, Qi an ge(x, %o-, t) satisfies a backward Fokker- 
Planck equation [201 



d t Q\ a nge{x, x , t)=D dlQi ange (x, x , t)-a D sgn(x - x ) d x Qi ange (x, x , t) (14) 

which holds for all x > with the initial condition Qi ang e(x, x ,t = 0) = 1 for all x > 
and the boundary conditions: (i) Qiange(^ = 0,Xo,t) = for all t (absorbing boundary 
at the origin) and (ii) Qiangc^ oo,Xo,t) = 1. 

To solve Eq. (14), we consider the Laplace transform Q\ an ge(x, xq, s) = 
Jo°° Qiange(a ; , ^o, t) e~ st dt, which, taking into account the initial condition, satisfies 

- 1 + s Q\ aJlge (x, x , s) = D - ^, a 2 ° ge _ a ° D sgn ( x ~ x °) ^dx^ ' ^ 

Making a shift Q\a, nge (x, %o,s) = 1/s + F(x, xo, s), one finds a homogeneous differential 
equation for F in x > 

d 2 F dF 
D -—^ - a -Dsgn(x - x ) — - sF = (16) 

with the boundary conditions: (i) F(x = 0, xq, s) = —1/s and (ii) F(x — > oo, xq, s) = 0. 
Taking the boundary condition (ii) into account, we obtain 

F(x, Xo, s) = Ai exp[— \i 2 (x — x )] forx>x (17) 
F(x, x , s) = Bi exp [— fii(x — x )} + B 2 exp [n 2 (x — x )] for x < x (18) 

where 



A*i = {y<Xo + 4s/D + ao)/2, and /x 2 = (Ja 2 + 4s/D - a )/2 . (19) 

From the boundary condition (i) at x = and using that F(x,Xq,s) and its first 
derivative d x F(x,x ,s) are continuous at x — x , the three unknown constants A%, 
B\ and B 2 are determined 

B x — 2//2 A x ; B^ ^-^ A,- andA^- 1 ^ + ^ 



(20) 

Finally, inserting x = xq leads to 



giange(a;0,s) = Qlangc^ = Xq, Xq, s) = - 



(fiX + fl 2 ) 



[2/j 2 e^ xo + (fii - fi 2 )e~^ x o] 
The mean first-passage time is then given by 

poo poo 

Ti ange (x , r) = t (-d t qi ange (x , t)) dt = qi angc (xo, t) dt = q\ ange (x , s = 0) . 
Jo Jo 



(21) 
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Taking the s — > limit in Eq. (21) leads us to 

1 



a 2 Q D 



[2(e 



a x 



1) - a x ] 



2(e 7 - r 



7 



r 



x 

D 



(22) 



where 7 = aox^ = \/r/D xq is the same dimensionless parameter as defined above (see 
Eq. Q). We can now compare the result in Eq. (22) with that of the resetting case, 
Eq. ([7]). Using the fact that e 7 — 1 > 7, we see that Ti ange (xo, r) > T reset (xo,r) for fixed 
values of x and D. 

The minimum value of Ti ange (xo, r) is obtained (see Fig. Q) for 7 = 72 where 72 
is given by the solution of <9 7 T lange (a;o, r) = 0, i.e., it is the positive root of the equation 
2(72—2) e 72 +72+4 = 0, leading to 72 = 1.24468 .... Thus the optimal mean first-passage 
time with Langevin dynamics of the searcher is given by 



rpOpt 

lange 



2(e 72 - 1) - 72 



D 



2.38762 



D 



(23) 



Comparing with the corresponding result in Eq. (|8j) for the reset dynamics, we conclude 
that 

1.54414... , , 



ry-iOpt 

^ reset 

ry-iOpt 

lange 



2.38762 



0.646728 



< 1 



This shows that the search process via the nonequilibrium diffusion combined with reset 
mechanism is significantly more efficient than the equilibrium Langevin dynamics of the 
searcher although the stationary distributions induced by both dynamics are the same: 
by this measure the nonequilibrium strategy beats the equilibrium dynamics. 



4. Multiparticle problem: Nonequilibrium vs. Langevin dynamics 

In the multiparticle version of the search process, we have a single immobile target at 
the origin and a team of searchers which are initially uniformly distributed on the line 
with uniform density p (see [23] for a colorful example). The searchers are independent 
of each other and the position of each searcher evolves stochastically (identical in law 
for all searchers) starting at its own initial position. This stochastic process, for the 
moment, is general. When any of the searchers finds the target, the search process is 
terminated. Let P s (t) denote the survival probability of the target up to time t, i.e, the 
probability that the target has not been found up to t by any of the searchers. To set the 
problem, we consider N searchers initially distributed uniformly in a box [— L/2, L/2] 
of size L and will eventually take the limit iV — > 00, L — > 00, but keeping the density 
p = N/L fixed. Let Xi denote the initial position of the z-th searcher. Thus, Xj is a 
random variable uniformly distributed in the box [—L/2, L/2]. Since the searchers are 
independent, it follows that for a given set of initial positions of the searchers {xi}, 



N 

Ps(t) = Y[q(xi,t) 
i=i 



(25) 
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Figure 2. Mean first-passage time T plotted versus the dimensionless parameter 7 
(setting xq = 1 and D = 1) for the nonequilibrium reset dynamics (shown by the 
solid (black) line) (see Eq. Q) and the equilibrium Langevin dynamics (dashed (red) 
line) (see Eq. ( |22[ )). Evidently, the optimal (minimum) T is higher in the equilibrium 
Langevin case. 



where q(xi,t) is the survival probability of the target due to a single searcher starting 
at Xj. For both (i) diffusion with resetting dynamics and (ii) Langevin dynamics in 
an effective potential V e s(x) = «o \x — xo\, the Laplace transforms of q(xi,t) are given 
respectively in Eqs. (J5]) and (21). The explicit dependence of P s (t) on Xj's has been 



suppressed in Eq. (25) for notational convenience. 



We first compute the average survival probability (averaged over the initial positions 



of the searchers) following Ref. [13J. Taking the average of Eq. (25) we obtain 



N 



N 

n 

i=l 



— > exp 



L 



L/2 



L/2 



(1 - q(xi,t)) dxi 



-P 



/oo 
i 1 
-oo 



q(x, £)] dx 



(26) 



where in the last step we have exponentiated for large L and then taken the 
thermodynamic limit iV — > 00, L — > 00 but keeping the ratio p = N/L fixed. Using 
further the symmetry q(x,t) = q(—x,t), one finally obtains the general expression of 
the average survival probability in the multiparticle system 



(P s {t)) = exp 



POO 

-2p [l-q(x,t)]dx 
Jo 



(27) 
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Actually, this result is rather general (see e.g. [19]) and holds for any stochastic process, 
the only assumption being that all searchers undergo the same stochastic process. If 
one can estimate the survival probability q(x, t) of the target at the origin for a single 
searcher starting at x, one has an exact formula for the average survival probability in 
the multiparticle problem. The quantity M(t) = J" °°[l — q(x, t)] dx has the interpretation 
as the expected maximum up to time t of the single particle process (starting at the 



origin) and the general formula in Eq. (27) and its discrete-time analogue have been used 
recently to compute exactly the average survival probability for a number of processes 
including continuous-time random walks as well as discrete-time Levy flights [23]. 

Similarly, the typical survival probability can be estimated [13] by first taking 
logarithm of Eq. (25), followed by averaging over the initial positions and finally 
reexponent iat ing 

N 



P^(t) 



exp 



1=1 



— > exp 



2p / \n[q(x,t)]dx 



(28) 



Thus, if we know the survival probability q(x,t) in the single searcher case, we can use 



the two exact formulae in Eqs. (27) and (28) to estimate respectively the average and 



the typical survival probability of the target in the multiparticle case. 
4-1. Average survival probability: nonequilibrium vs. equilibrium 



Here, we analyze the asymptotic large t behavior of (P s (t)) in Eq. (27) and compare the 
expression obtained for the resetting dynamics with the one for the Langevin dynamics. 
Let us denote 



M(t) 



[l-q(x,t))dx so that (P s (t)) = exp [-2pM(t)\ 



(29) 



In terms of the Laplace transform of M(t), M(s) = J °° M(t) e st dt, Eq. (29) reads 

poo |~ -| 

M{s) 



POO 


- ~ q(x,s) 


Jo 


s 



dx 



(30) 



We now consider the two cases (i) diffusion with reset and (ii) Langevin dynamics 
separately. 



4-1.1. Nonequilibrium case (diffusion with reset) Substituting q rese t(x,s) from Eq. (JSJ) 

(31) 



in Eq. (30) and performing the integration over x, one obtains exactly [13 

M(s. 



jDir + s) , / V + s 
— In 1 



r s V s 

The large t behavior of M(t) will be derived from small s behavior of M(s): 

Id 



As s -> . 




Ins ln(r) 



(32) 
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Inverting the Laplace transform, we find that the leading asymptotic behavior of M(t) 
for large t is given by 



M(t) 



hit + \ — ln(r) + . . . 



(33) 



Thus, the average survival probability from Eq. (29) decays algebraically for large t 



(P r s e (t)) ~ oi r" 1 where 1 = 2p v / D/r and en = exp {-6 X lnr 



(34) 



4-1-2. Equilibrium case (Langevin dynamics) Here, we use the expression (21) of 



Qungeix, s) in Eq. (30), perform the integration over x exactly and get 



M{s) 



Hi + Hi 



1 



/'i 



2/ii + fl 2 fi2 - Hi 



(35) 



_2p 1 jj, 2 s\ V H1+H2 Hi + A*2 2/i 2 
where /ii j2 are defined in Eq. (19), a in Eq. ([3]) and the hypergeometric function 
2-Fi (a, 6, c, z) is given by 

, N a6 a(a + l)b(b + l) z 2 a(a + l)(a + 2)6(6 + 1)(6 + 2) z 3 

2 Fi(a, b,c,z) = l-\ z + — ' \ — + — —, — + . . . 

u ' ' ' J c c(c+l) 2! c(c+l)(c + 2) 3! 

(36) 

Expanding /i 12 for small s, keeping terms up to 0(s) and using a^D = r, one gets 



M(s) 



cto D 



, 2s 
1 + — + 

r 



2*1 



1,1-^,2-^,-f -i + A 
r r 2s 2 2r 



(37) 



To make further progress, we use the following identity 2^(1,1,2,-^) = ln(l + z)/z 



which follows from the definition Eq. (|36|). Expanding for small s, we find the following 

ln(r/2) + . . . 



leading order behavior of M(s) from Eq. (37) 

M(s) -- 



1 i / \ 1 
— In s + — 



(3* 



aos ao s 

where . . . correspond to lower order terms that vanish as s — > 0. This indicates that for 
large t, 

M(t) w — ln(t) + — ln(r/2) + . . . (39) 



Hence, from Eq. (29), we obtain the average survival probability 

(P] angc (t)) w a 2 r 82 where 6 2 = 6 1 = 2 p ^fD~jr~ and a 2 = e~ dl ln{r/2) 



(r/2) 



-0i 



(40) 

Thus the power law exponent 9 2 characterizing the algebraic decay of the average 
survival probability in the Langevin case is identical to the nonequilibrium case, though 
the amplitude a 2 = 2 01 a\ is larger than a±. Hence, for large t, the average survival 
probability in the Langevin case is greater than that of the nonequilibrium case 



(preset ^ < (P^^)) . 



(41) 



Thus, we conclude that on average the target is found faster in the nonequilibrium case 
than in the equilibrium one. 
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4-2. Typical survival probability: nonequilibrium vs. equilibrium 

We now analyze the asymptotic large t behavior of the typical survival probability, 



defined in Eq. (28). We define 



POO 

W(t) = / ln[q(x, t)} dx so that P s typ (t) = exp [2 p W{t)\ . 
Jo 



(42) 



We recall that the Laplace transform of q(x, t) is denoted by q(x, s) and its expressions 
are given in Eqs. ^ and (21) respectively for (i) the diffusion with reset dynamics and 
(ii) the equilibrium Langevin dynamics in the effective potential V e s(x) = ocq \x — xq\. 
The two cases will be considered separately. 

4-2.1. Nonequilibrium case (diffusion with reset) To analyze the asymptotic large t 
behavior of W(t), we need to know how q(x, t) behaves for large t. The Laplace transform 
of q(x, t), given in Eq. (|5]), has a pole at s — (for fixed r and fixed x) which satisfies 



Si + r exp 



r + Si 



D 







(43) 



Clearly s± = si(x) depends implicitly on x and one must have si(x) < 0. Then, to 
leading order for large t, it follows from Laplace inversion that 



q(x, t) ~ exp [si(x) t] = exp [— |si(x)| t] . 



Consequently, from Eq. ( 42 ) , we find to leading order for large t 

' POO 

W(t)pa- / |si(a;)|dac 
Jo 

Hence, the typical survival probability decays exponentially for large t as 

POO 

Pl yp (t) ~ exp [— 2 p Ki t] where /ti = / \si(x)\dx. 

Jo 



(44) 



(45) 



(46) 



To compute si(x), it is useful to first define si(x) = — r(l — u), so that Eq. (43) 
reads, in terms of u and the dimensionless length z = «o x 



u — 1-|- exp [— y/uz\ = . 
Clearly, as z — > oo, u(z) — > 1 and as z — > 0, u(z) —> 0. Hence 



|si(x)| dx 



[1 — u(z)] dz . 



(47) 



(4J 



'o a o Jo 

The idea then is to transform the integral over z to an integral over u. We then use 
dz = du/\u'(z)\ where u'(z) = du(z)/dz. The derivative can be easily computed from 



Eq. (|47j). Expressing Eq. (|47j) as z = — ln(l — u)/ y/u and taking derivative, we get 

1 



dz 



du 2u*l 2 



ln(l - u) 



1 



Mil — u) 



(49) 
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Hence this gives, using a = y/r/D, 



poo p 1 

Kl = y/r~D / [l-u(z)]dz =VrD / (1 
Jo Jo 



u 



2u 3 / 2 



ln(l -u) + 



u (1 — u) 



The integral in Eq. (50) can be done explicitly to give 



«i = 4(l-ln2)Vr\D. 

Hence finally we find 

Pf p {t) ~ exp [-8 (1 - \n2) pVrDt 



cxp 



- (2.45482... )pVrDt 



du . 
(50) 

(51) 
(52) 



As explained in [13] . the fact that the average and the typical survival probabilities 
have different behaviours in the large time limit is a consequence of the memory of the 
initial conditions in the diffusion process with resetting. 

4-2.2. Equilibrium case (Langevin dynamics) In this case, we proceed in exactly the 
same way as the nonequilibrium case, except that to evaluate W(t) = J °° ln[g(x,t)] dx 



in Eq. (42), we need the expression (21) of the Laplace transform of q\ ange (x, t). The 
function %,nge(^ ', s) has a pole at s — S2(x) which is a root of the equation (with fixed 
r and fixed x) 



2p 2 {s 2 ) exp[//i(s 2 )x] + (/ii(s 2 ) -/U 2 (s 2 )) exp[-// 2 (s 2 ) x] =0 



(53) 



where we emphasize that Pi^(s) defined in Eq. (19) depend on s. Because /ii > /i 2 



one must have /i 2 (s 2 ) < and therefore s 2 (x) < to fulfil this equation. Thus, we find 
q(x,t) ~ exp[— |s 2 (x)| t] for large t and W(t) ~ — [J °° |s 2 (x)| dx\ t. Consequently 



where 



P* yp (t) = exp [2 p W{t)] ~ exp [-2 p k 2 t] 

POO 

K>2= I \S2(X) \ dx . 



(54) 



(55) 



To compute k 2 , we reorganize Eq. (53) slightly (using the explicit expressions of p\ and 



/i 2 ) and express it in terms of the dimensionless length z = a^x as 

0. 



1 + 4s 2 /r — 1 + exp — A/l + 4s 2 /r2 



(56) 



where we have used a^D = r. Let us further define v — 1 + 4s 2 /r in terms of which Eq. 



(56) reads 



y/v — 1 + exp [— y^v z] = . 



Hence, from Eq. (55) we get 



«2 



[1 — v(z)] dz 



(57) 



(58) 



Optimal diffusive search: nonequilibrium resetting versus equilibrium dynamics 14 



where v(z) is the solution of Eq. (57). Now the solution of Eq. (57) has two branches: 
v = 0, < z < 1 and a branch with v(z) > for z > 1. Thus (58) becomes 



«2 



f(z)] dz 



(59) 



To compute the integral in (|59|), we use the same trick as above. 

dz(v) 



/oo r\ 
[i - v(z)\ dz = [i 



-v 



dv 



dv 



(60) 



where z (v) = — ln(l — \/v)/y/v, which follows from (57). Here, the integral (60) is easily 
evaluated by integration by parts 

• i 



\l-v)^-dv= [{l-v)z(v)t =0 + f z(v)dv 



- 1 - 



1 ln(l 



,1/2 



dv — 1 . 



Therefore from (59) we have 



Hence finally we obtain 



K 2 



pty P{t) 



exp 



-p v r.D t 



(61) 



(62) 



(63) 



Comparing this result to Eq. (52), we find that for the multiparticle well, 
the typical survival probability of the target in the equilibrium Langevin case is larger 
than the nonequilibrium dynamics of diffusion with resetting 



langc 



(64) 



This means that the target is found faster in the nonequilibrium case than the 
equilibrium Langevin case. In other words, the target search process by multiple 
searchers, as in the case of a single searcher, is more efficient when the searcher undergoes 
nonequilibrium diffusion and reset dynamics, rather than the equilibrium Langevin 
dynamics. 



5. Concluding Remarks 

Equilibrium thermodynamics teaches us that dissipation is reduced when a system 
remains close to an equilibrium state and that the transformations that affect it are 
quasistatic and reversible. Such a statement of local optimum can not be taken as 
a paradigm: it can be advantageous in some circumstances to be driven away from 
equilibrium, by creating non-vanishing stationary currents that break detailed balance 
and time-reversal invariance. 

Optimal search problems provide us with concrete examples in which nonequi- 
librium can defeat equilibrium. In the present work, we have undertaken a systematic 
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comparison between two related search strategies. We consider a predator who performs 
an intermittent search by alternating phases of free diffusion and resetting jumps. The 
invariant distribution corresponding to such a dynamics is a nonequilibrium stationary 
state in which a non- vanishing probability current is directed towards the resetting posi- 
tion. However, one can readily define a fictitious equilibrium Langevin dynamics leading 
to the same stationary distribution (but with different local transition rates, leading to a 
vanishing current in the steady state). For the resetting model, this amounts to defining 
an effective attractive potential centered on the resetting position. The hitting-times 
to an immobile target can be exactly calculated in both dynamics (nonequilibrium vs 
equilibrium) and compared. The same problem can be addressed in the case of a team 
of independent searchers; here the typical and the average survival probabilities differ. 
Nevertheless, in all the cases we have studied in one dimension, we have shown that the 
nonequilibrium strategy is advantageous: the intermittent search is always more efficient 
than equilibrium diffusion. It would be interesting to see if the same conclusion holds 
in higher dimensions as well and in more realistic search scenarios where, for example, 
there is a time penalty for resetting. 
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